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Abstract 

Transport through generalized trees is considered. Trees contain the simple nodes and supernodes, either 
well-structured regular subgraphs or those with many triangles. We observe a superdiffusion for the highly 
connected nodes while it is Brownian for the rest of the nodes. Transport within a supernode is affected by 
the finite size effects vanishing as A'^ ^ oo. For the even dimensions of space, d = 2, 4, 6, . . ., the finite size 
effects break down the perturbation theory at small scales and can be regularized by using the heat-kernel 
expansion. 

PACS codes: 87.15.Vv, 89.75.Hc, 89.75.Fb 

Keywords: Structures and organization in complex systems, Networks and genealogical trees. Diffusion 

1 Introduction 

Complex networks have become a new paradigm in physics. They have been studied exten- 
sively due to their relevance to many real systems from the World Wide Web to the biological 
and social networks. Graph theory has been successfully applied to a wide range of different 
disciplines requiring a description of sets of elements either connected or interacting pairwise. 
Geometry and topology of underlying graphs have a deep influence on the physical properties 
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of complex networks. In view of that one is interested in the properties of graphs which affect 
the dynamical behavior in the models defined on them. On the other hand, the study of large 
complex systems calls for the statistical methods, to give an effective description of collective 
dynamical behavior observed in them. 

In the present paper, we study the nonlinear transport on large inhomogeneous graphs. 

Diffusion processes defined on the 1-dimensional small world networks have been considered 
in [Ij. It appears that the mean field theory breaks down for the small worlds models in 
dimensions (i < 2 [2J due to the emergence of strong site-to-site fluctuations of the Green's 
functions. In [3], the e-expansion for a small world model has been developed where e = 
2 — d (quantifies the departure of the space dimension from 2). The small world network is 
constructed by adding random shortcuts to a regular lattice, and the density of links p > has 
been taken enough small: pa'^ <^ 1 where a is the lattice scale. The Green's functions have 
been computed for d = 2, but the e-expansion breaks down for e = 1 due to traps appearing 
in the system at d = 1. In contrast to the small worlds model [3], we consider the diffusion 
processes defined on graphs approaching the networks known in the literature as "froths" [1]; 
they could contain the regular subgraphs but are, in general, far from being the regular lattices. 
The space tiled by the froth can be curved, and, in such a case, the intrinsic dimension of the 
cellular system does not coincide with the dimension of the embedding space. 

Transport phenomena had been extensively studied for the disordered media [S]-IZ] and 
fractals [8]- [10]. It has been revealed that the transport properties depend upon the intrinsic 
dimension df of disordered systems. In particular, the random walk probability P is propor- 
tional to t~^f ^'^ f {x / ^/t) where / is a scaling function. Similar results have been reported for 
nonlinear diffusion processes |llj-|16j. 

In our model, we consider the complex networks as the generalized trees with two types of 
nodes: the simple nodes (more often of low connectivity) and the supernodes which are either 
the subgraphs containing many polygons or the fc-regular subgraphs. In a continuous setting, 
the supernodes can be treated as complete compact curved Riemann surfaces characterized by 
the finite areas (see Sec. 3). The supernodes are bridged by the tree components. We study the 
large scale asymptotic behavior of nonlinear transport process defined on such the generalized 
trees. It is affected by both the space curvature within a supernode and the variation of space 
dimension between the supernodes. 

Strictly speaking, the possibility to replace a finite graph by a compact continuous manifold 
for a nonlinear diffusion process is a challenging question. For the linear differential operators 
acting on the periodic functions, it requires that the spectra associated to the discrete and 
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continuous problems are similar [T7] . In particular, the multiplicities of the eigenvalues should 
be equal for both problems. 

However, by this time, it is still not much known even about the spectral properties of 
Schrodinger operators defined on curved compact manifolds. Some estimations on the multi- 
plicities of eigenvalues and the Euler characteristic of a surface can be found for the 2D-sphere, 
for the Klein's bottle, and for the 2D-torus in [l8]-[22]. Therefore, even in this case, the justi- 
fication of the validity for the approximation of processes defined on the graphs by that ones 
defined on the compact continuous manifolds is a difficult problem. Concerning the nonlinear 
diffusion process, we should confess that such an approximation is always an assumption. 

We consider the diffusion as a generalized Brownian motion with arbitrary boundary condi- 
tions described by a functional integral (see Sec. 4). The nonlinear term included into the diffu- 
sion equation models the effect of varying space dimension, the possible fluctuations of transport 
coefficient, the diffusion-reaction processes, and the possible queueing due to a bounded trans- 
port capacity of edges bridging the supernodes (see Sec. 5). The fluctuations have been treated 
in the framework of fleld theory approach. 

The large scale asymptotic behavior of Green function is a superdiffusion for the nodes 
of high connectivity, while it is still Brownian for the rim nodes (see Sec. 6). In the regular 
subgraphs endowed with the standard orientation of edges, the effect of space " curvature" within 
a supernode rises the flnite size corrections to the scaling behavior (see Sec. 6 and Appendix 
A). In Sec. 7, we use the heat-kernel expansion to regularize the effective action functional at 
the small scales for the even dimensions of space (i = 2,4,6,.... 



2 Description of the model and the results 

We are interested in the large scale asymptotic behavior of Green function for the diffusion 
equation for the density u{x,t) deflned on the curved Riemannian manifold with metric tensor 
Qij, with the nonlinear term oc (a > 1) included. We assume the such a model can be 
considered as an approximation of nonlinear diffusion process deflned on the graphs containing 
the large regular subgraphs. This approximation, however, is rather intuitive and cannot be 
proven nowadays rigorously. If the assumption is true, then logg deg(x) plays the role of effective 
local space dimension at node x on the graph, deg(a;) being the degree of x. The nonlinear 
diffusion term is relevant to the large scale asymptotic behavior if 

2 

log2 deg(x) 

and is irrelevant otherwise. 
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In the critical phenomena theory, the physical degrees of freedom are replaced by the scaling 
ones related to each other by the RG transformations. The large scale asymptotic behavior of 
physical system is then determined by the properties of scaling invariant system at a certain 
stable fixed point of the RG differential equation. The derived results are irrelevant to the order 
of space and time (infinite) limits. 

Our main result is that for the large scale asymptote of Green function, 



The main technical difficulty of computations in the curved space metrics is the lack of the 
global frequency-momentum representation. In order to regularize the effective action at small 
scales in the one-loop order, we have used the heat-kernel expansion which does not depend 
upon the space-time topology. In the regular graphs endowed with the standard orientation of 
edges (the cycling ordering of edges to be the same for all nodes), the effect of curvature within 
the supernodes reveals itself by the finite size corrections to the scaling behavior. 

3 The regular subgraphs viewed as Riemann surfaces 

In the present paper, we consider the fiat Euclidean space R*^ as the limit a ^ of a regular 
lattice La = with lattice scale length a. For simulating the diffusion equation ut = Au for 
the scalar function u defined on the lattice, one uses its discrete representation. 



where Ux is the lattice neighborhood of a; G The cardinal number 2'^ is uniform for a given 
Ca, and d is interpreted as the dimension of Euclidean space. 

Being defined on an arbitrary connected graph F = (Vr,-Er)5 in which Vr is the set of its 
vertices and Er is the set of edges linking them, the discrete Laplace operator has actually the 
same form as in ([1]) excepting for the cardinality number changed to 2^=^ where 6x = log2 deg(a;) 
and deg(x) is the degree of vertex x G F. By counting the number of vertices |Vr|, the number 
of edges \Er\, and the number of faces |Fr|, one obtains the Euler characteristic of a planar 
graph F, 





(1) 



X(F) = |Vr| - l^rl + l^r 



and its genus. 



g(r) 



2 - x(r) 



(2) 
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the number of nonintersecting cycles on the graph T. We are interested in the subgraphs of T 
ample with polygons, in particular, the triangles, which can be used as a base for a local mesh 
of a Riemann surface (the Delaunay triangulation). 

An intuition on the deep relation between the triangle structures in the complex networks and 
their curvature has been expressed in "Triangles capture transitivity, which we measure by 
the associated notion of curvature". In [27], upper and lower bounds on the number of graphs 
of fixed degree which have a positive density of triangles have been estimated. In particular, it 
has been shown that the triangles seem to cluster even at low density. Then, in [28] it has been 
proven that the probability for a randomly selected vertex to participate in T triangles decays 
with T following a power-law, if the graph is scale free with the degree exponent satisfying 
2 < (3 < 2.5. Moreover, if /? = 2 + 1/3, the density of triangles appears to be finite in such 
graphs. 

The Laplace-Beltrami operator on a triangle mesh has been defined in connection with the 
various graphical applications such as mesh fairing, smoothing, surface editing in 3Z)-space 
[29]-[32]. The stiffness matrix correspondent to it can be computed for each triangle specified 
by its vertices x = (xi, X2, X3), y = {yi, 1/2,113), z = {zi, Z2, Z3) as 

f X2-X3, y2- Z2- Z3\ 

X3 -xi 1/3 - yi ^3 - Zl 



K = 5aBB^, B = -i 



A 



\ xi- X2 yi- y2 Zl- Z2 J 



where S'a is the area of triangle. 

However, from the various models of complex networks, it seems rather easier to count 
the number of edges which link a node to others than to check out if its neighbors are really 
connected forming triangles. Instead of hunting for triangles, while analyzing the graphs of real 
world networks, one can search for the 3-regular subgraphs srF for which there are precisely 
three edges incident at each vertex (loops are counted twice, multiple edges are allowed). These 
graphs form a honeycomb. The idea of using 3-regular subgraphs to study the topological 
properties of rendered Riemann surfaces has been proposed in |33j. In [Ml ES], it has been 
shown that for each 3-regular graph s^jF of 2N nodes with an orientation O (which assigns 
to each vertex of a/jF a cyclic ordering of edges incident at it) one can construct a complete 
Riemann surface S^^rT, O) by associating the ideal hyperbolic triangles to each vertex of a/jF 
and gluing sides together according to the edges of the graph ^rV and the orientation O. 
The resulting surface is endowed with a metric excepting for the finitely many points (cusps) 
where the metric could be undefined. The surface area is finite and equals to 2t{N . The 
conformal compactification of S{3j{T, O) are dense in the space of Riemann surfaces. The 
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standard orientation on the 3-regular graph which is the 1-skeleton of a cube contains six left- 
hand-turn paths on which a traveller always turns left, giving that the associated surface is a 
sphere with six punctures. The choice of different orientations can give surfaces of genus 0,1, 
and 2, [31]. 

The closed paths of length k on are homotopic to the closed geodesic lines on S'(r, O). 
Their cardinality is defined by the spectral density of aRF-subgraph that is the density of 
eigenvalues of its adjacency matrix [23] . 

which converges to a continuous function with N —>■ oo [Xj is the j— th largest eigenvalue of 
adjacency matrix A). The k—th moment Mk of g{X) scales with the number of cycles Ck of 
length k, 

Mfc = - y = -Tr(A^) = -Ck, 



Ck can be computed by using of the zeta-function, 

^^^^ " det(l-;2A)' 

Namely, 

1 d'^ 
^'=Jk^\dz^^^'\^, 
All roots of the polynomial det(l — zA) lay on the unit circle and coincide with the inverse 
of the eigenvalues Aj^ of adjacency matrix excepting for the zeros, {A = 0}. It is worth to 
mention that S{^^ , O) allows for a conformal compactification. 



4 Model of nonlinear diffusion through complex networks 

In the previous section, we have supposed that a regular subgraph in a complex network can 
be treated, in a continuous setting, as a Riemann space of finite area characterized by some 
metric tensor. For a particular time slice, the line element ds between each pair of neighboring 
points on the spatial surface is given by 

ds^ = gij{x)dx^dx^ 

ij 

where dx'' denote the differences between neighboring points, and Qij denotes the metric tensor. 
Then, the complex network as a whole can be considered as a disordered media in which the 



6 



compact islands S of curved Riemann space (with the "effective" space dimensions Ss) are 
bridged by the tree hke graph components (see Fig. 1) in which the local space dimension 6x 
can vary from point to point. 

The transport properties through such a disordered media is essentially of nonlinear nature. 
In the previous studies of nonlinear diffusion [TT] -[15], the authors had introduced various 
nonlinear terms into the diffusion equation modelling the possible fluctuations of transport co- 
efficient, the diffusion-reaction processes, and the queuing due to a bounded transport capacity 
of edges. We also introduce it for accounting the effect of varying dimension of space in the 
complex network (see the discussion below). To be certain, let us consider the equation for the 
scalar density field w(x, t) defined on a Riemann surface, 

d'^u d 
Vtu = 9o^O g^ - ^o^oRu - Vo^ou", '^* = '^* + ^o^- (3) 

All 0-subscripted variables denote their bare values before the application of the renormalization 
group transformation. In Eq. ([3]), the Riemann metric tensor depends upon the chosen 
conformal parameterization of regular subgraphs, R is the scalar curvature. We prefer to 
keep the entries of Qq dimensionless therefore we have introduced the parameter z/q having 
the dimension of a viscosity. For the purpose of this paper, we will examine metric rescalings 
that are spacetime constants (we suppose that the subgraph of F is regular and the edges do 
not rewire with time). However, it is possible to consider the effect of a rescaling given by a 
spacetime-dependent function. 

The covariant derivative Vt contains the curvature drift term proportional to 6q = 5'o'^F*^ (the 
curvature drift velocity) which expresses the local anisotropy of space because of its curvature. 
F*;, are the Christoffel symbols calculated out from the metric tensor g'^^ in the standard way, 
— d^Kdkij + 9ji,k — gjk,i)/'^- While being interested in the long time large scale ranges, one 
usually keeps only the first order derivative term di if it presents in the linear part of equation 
since its contribution 0{k) should dominate over the diffusion 0{k'^) for small k. Nevertheless, 
we keep both terms to ensure the convergence of integrals in time, in the limit of flat metric. 
However, in the general curved spacetime, the homogeneity required for the existence of a global 
momentum space representation is lacking. 

The real valued parameters and rjo are the bare coupling constants governing the coupling 
of configurations u to the scalar curvature of space R and to the varying effective space dimen- 
sionality Sx respectively. We assume them to be small. The nonlinearity exponent is a > 1 (a 
is not necessary integer) [11] -[16]. 

Let us explain the role played by the nonlinear term in ([3]) in more details. In the critical 
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phenomena theory, the physical degrees of freedom are replaced by the scaling degrees of free- 
dom. In particular, one considers the canonical dimension dp instead of time and space physical 
dimensions of the quantity F. The analysis of canonical dimensions allows for the selection of 
relevant interactions among all possible interaction which could arise in the model. In the spirit 
of critical phenomena theory, the nonlinear term that could affects substantially the large scale 
asymptotic behavior typical for the diffusion process should have the same canonical dimension 
as the normal diffusion. If the canonical dimension of the nonlinear term added to the equation 
is less than of the ordinary diffusion, it has to be neglected. The opposite is also true: if the 
nonlinear term provides the leading contribution to the asymptotic behavior, then the diffusion 
term has to be dropped. 

Therefore, it is interesting to consider the model in which the nonlinear term would play 
an important role. Below, we demonstrate that the exponent a is related to the dimension of 
space d and the dimension [r^o] of coupling constant r/o- If we assume for a moment that in the 
vicinity of some point the dimension of space d is changed to some other value 5, then, strictly 
speaking, Eq. (3) could have no sense therein: either the diffusion term or the nonlinearity 
should be neglected. For given a and [?7o], there is only one value d at which the Eq. (3) 
is relevant with respect to the large scale asymptotic behavior of diffusion process. If we 
consider the plane of parameters [r/o] and a, then the relevant space dimensions d is a line on 
it. Therefore, by tuning the values of a and [r^o], one can "modify" the space dimension d in 
the model of nonlinear diffusion. It is indeed unphysical to change the nonlinearity exponent a 
(we suppose that a is a property of a certain physical process), however, one can tune the value 
of [r^o] and use it as the small expansion parameter of perturbation theory (like the parameter 
e = A — d in the Wilson's theory of critical phenomena). 

A similar idea is used in the usual dimensional regularization of Feynman diagrams. In the 
continuous Euclidean space, the dimensional regularization scheme does not look natural and 
therefore is usually treated as a formal trick which helps to reformulate the singularities arisen 
in the Feynman graphs in the form of poles in e. However, if the dimension of physical space 
could vary, than the dimensional regularization would acquire the natural meaning provided the 
nonlinearity exponent a is fixed and the correspondent nonlinear diffusion process is relevant 
with respect to the large scale asymptotic behavior. 

As we have mentioned above, such a relevance can be justified by means of dimensional 
analysis [TB] . Dynamical models have two scales: the length scale L and the time scale 
T. The physical dimension of viscosity is [z/q] = L'^T~^, of the scalar curvature is [R] = L~^, 
and of the drift velocity is [6q] = LT~^ . Let us choose the physical dimension for the coupling 
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constant to be [?7o] = L"^'^ (assuming that e = in the logarithmic theory when the nonhnear 
interaction is marginal) and note that the dimension of a scalar quantity is [u] = L~'^. In the 
free theory, dt oc and in the analysis of canonical dimensions one has ~ T. All terms 
in ([3]) should be of the same canonical dimension, in particular [dtu] = [tiqUqu"'], therefore 
—2 + logj;^['u] = Q;log^[M] — 26 and finally, [u] = L'^'y^^'^)l'y^-°'\ from which it follows that 

d = 2.\^. (4) 

1 — a 

The above relation gives us a hint of that the space dimension in the model of a nonlinear 
diffusion can be effectively tuned by the parameters a and e. We choose the parameter e to 
quantify the local irregularity of the graph by measuring the relative deviation of the node 
degree 2^"^ from the cardinality number 2^^ in the regular lattice, 

.= l-f (5) 

Ox 

The nodes with deg(x) < 2'^ correspond to £ < 0, while the nodes for which deg(x) > 2'^ are 
described by £ > 0. 

We supply Eq. ([3]) with the locally integrable initial condition m(x, 0) and study the standard 
Cauchy problem being interested in the large scale asymptotic Green's functions ^(x, xq; t, to)- 
For a curved space, the natural way to proceed is to examine the change to the Green's functions 
as the metric is scaled. This can be achieved by moving the points along the geodesies (cycles 
in the graph F) connecting them or alternatively by scaling the geodesic distance function (the 
metric). 

The Green's functions of nonlinear problem ([3]) supplied with the integrable initial conditions 
can be formally calculated by the perturbation series with respect to the nonlinearity (as the 
coupling parameter rjQ is small) followed by the integrations over the initial condition ^(x, 0). 
Some integrals estimating corrections to the linearized diffusion problem diverge logarithmically 
since the integration domain is not compact. If we introduce the e-parameter in accordance to 
dS]), the divergences reveal themselves by the poles in 

Therefore, the nonlinear interaction is irrelevant (in the sense of Wilson) for e < (rf > 
2/ {a — 1)), but is essential as e > when the ordinary perturbation expansion (in the form of 
series in t^q) fails to give the correct large scale asymptotic behavior and the whole series has 
to be summed up. For instance, it happens at d = 2 for a = 2. 

In other words, the logarithmic (marginal) value of a is determined by comparison of the 
nonlinear contribution with that of the linear dissipative term, a log = 1 + 2/d. While introduc- 
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ing the parameter e accounting for the local change of connectivity in the graph, we effectively 
pass from a log to its new value a [^^ = 1 + 2/5x- Then it turns that the nonlinear contribution 
to the long range asymptotic transport through the rims {e < 0) is irrelevant in comparison 
with the linear diffusion and therefore can be neglected. In contrast to it, the contribution com- 
ing from hubs {e > 0) is more essential than the linear one and has to be taken into account in 
all orders of perturbation theory since the relevant fluctuations dominate the diffusion at large 
scales. 

We calculate the asymptotic Green's functions for the model ([3]) in the logarithmic theory 
(on the regular subgraphs with the cardinality number 2'^) in curved space metric and develop 
the e— expansion accounting for the corrections in the long time large scale region due to the 
irregularity of graph. The relevant contributions to the nonlinear transport coming from hubs 
are summed by the field-theoretic renormalization group method. Herewith, the real values of 
parameter e has to be taken as = 1 — d/6x, the excess of hub's connectivity over the regular 
cardinality number 2*^. 

We conclude this section by remarking that the problem of renormalization in a curved 
spacetime has been discussed extensively in the literature ([37]- [15] and by other authors). 
However, it has never been studied in connection with the critical phenomena theory. The 
analysis of transport through the graphs would provide us with such a model. It is worth to 
mention that in contrast to the case of the gravitational field, we are not restricted on graphs 
by the equivalence principle, so that the curved space has not to be flattened in any sufficiently 
small region. 



5 Diffusion as a generalized Brownian motion 

We set up a field theory formalism for the study of asymptotic properties of the Green's 
function by means of the renormalization group (RG) equations which are valid in the curved 
spaces. Our method is closely related to the approach discussed in [38] regarding the Green's 
functions as the functions of metric and scaling the metric instead of scaling the coordinates or 
moments. 

It is well known that many problems of stochastic dynamics (and of the transport through 
a disordered media, in particular) can be treated as a generalized Brownian motion, P{u) = 
{6 {u — u (x, t))), in which the average is taken over all configurations of field m(x, t) satisfying 
the dynamical equation 

Vtu - uq^lbu + i'oT]qu" + ioVoRu = -^S{t - to)5(x - Xo), (6) 
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for the integrable initial condition ^(x, 0) and arbitrary conditions on the boundaries of the 
graph. The Laplace-Beltrami operator Alb is given by ([3]), and g = \ detg^l- 

We use the functional representation of the 5— function for expressing the probability, 

in which u marks the position of a "particle", and the auxiliary field u' (of the same nature as 
u) is not inherent to the original model, but appears since we treat the dynamics as a Brownian 
motion. The formal convergence requires the field u to be real and the field u' to be purely 
imaginary. Should a unique solution of dynamic equation exists, we perform the natural change 
of variables in ([7]), 

{u - M(x,t)) -ViM + UqAlbU - UoT]ou" - ^qVoRu - to)^(x - Xo) = 0, 

V 9 



from which it follows that 



P{u) = j j 'Du'expS{u,u') det M 



S = So- T]oiyoTTg{u'u°') —u'{xo, h), = Tig {-u'VtU + uqu'Albu - ^qPqRu'u) . (8) 



where det M is the Jacobian associated to the change of variable, and S{u,u'), the action 
functional, 

1 

79' 

The trace Tr^ means the summation over the discrete indices and the integration J dvx J dt 
over the invariant volume element on the d-dimensional manifold, dvx = \/g{x)d'^x. 

The Jacobian det M deserves a thorough consideration. The linear part of the variable 
transformation can be factorized from it, 

det M = det Mq det ( 1 - A„„/ Mi ) (9) 

where Mq = —Vt + i^qAlb — ^o^oR, the interaction part Mi = az/o?7ou"~^(5(t — t'), and A^u' is 
the Feynman propagator in curved Riemann space, [1^ W7\ , defined as the solution of linearized 
problem 

(-Vt + z/qAlb - Coi^oR)^uu'{^ - x, t - t') = —5{t - t')6{^ - x'). 

V9 

For more details as well as its explicit form see Appendix A. 

It is important for us that the propagator A^u' is proportional to the Heaviside function 
9(t — f) as a consequence of causality principle. The first factor in ([9]) does not depend upon 
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fields and therefore can be scaled out of the functional integration. The second factor in 
can be expanded into the "diagram" series, 

logdet(l - A„„,Mi) = -Trg(A„„,Mi + ^A„„,MiA„„,Mi + . . .) 

comprising of cycles of retarded Feynman propagators proportional to the Heaviside functions 
and therefore being trivial, excepting for the very first term, Trc,(A„„'Mi), in which the operator 
Ml contains 6{t — t'). The first term in the expansion is proportional to the undefined quantity 
6{0) which value is usually taken as 1/2. However, in the critical dynamics, another convention 
is used [H], ^(0) = , under which the Jacobian det M turns to be just a constant and therefore 
may be scaled out away at the irrelevant cost of changing only the normalization. 

The action functional of type ([8]) in the problem of nonlinear diffusion in the fiat metric has 
been introduced in [16j. In [IS], the functional with an ultra-local interaction term like in ([H]) 
has been derived as a limiting one in the framework of MSR formalism (stochastic quantization, 
[50]). The renormalization of field theoretic models with ultralocal terms, located on surfaces, 
had been studied in [5T1| in details. 

Further insight into the field theory representation of Brownian motion and the properties 
of auxiliary field u' can be obtained from the equations for the saddle-point configurations. The 
first equation, 6S/6u' = 0, recovers the original Cauchy problem. The other one, 6S/6u = 0, 
reads as following 

Vtu' + uqAlbu' - ioVoR = arjououu" 

and is characterized by a negative viscosity. One can conclude from it that the auxiliary field 
should be trivial for positive time, u'{t > 0) = 0, and decays as t ^ —oo. 

In the framework of field theory approach, the Green function G{x,t; XQ,to) for the Cauchy 
problem ([6]) can be computed as the functional average, 

fVufVu'u{x,t)expS{u,u') 
G{x, t; xo, to) = {u) = j^,j^^.,,^^So) ' ^''^ 



in which S is the action functional ([8]). The Green function (iTOl) and all higher moments of fields 
u and u' allow for the standard Feynman diagram series expansions. The diagram technique 
with the ultralocal interaction terms has been discussed in [TB]. A special feature of such 
diagrams is that the final point of any diagram corresponds to (xo,to)- It is worth to mention 
that diagrams could formally contain a non-integer number of lines (since the nonlinearity 
exponent could deviate from an integer number). Diagrams are drawn of three elements: i) the 
final point (xo,to) with an arbitrary number of attached w'-legs (we mark them by a slash, see 
Fig. 2); ii) the interaction vertex with one -u'-leg and a M-legs attached to it (we put the letter 
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a inside the loop to stress that it is not necessarily integer); iii) the propagator A„„/ is only 
available. The first three diagrams for the Green function (fTOl) are sketched in Fig. 2. 

6 Renormalization Group equation and scaling behavior of scalar 
field coupled to a complex network 

Power-counting arguments [16] show that the divergent Feynman diagrams are those which 
involve any number of external w'-legs (lines with slashes). For all these functions the formal 
index of divergence equals to zero, therefore all divergent contributions are logarithmic (the 
correspondent counterterms are constants). It is worth to mention that the model is renormal- 
izable despite the fact that it requires an infinite number of counterterms. It is sufficient to 
renormalize the only "one-particle-irreducible" Green function (the only diagram block which 
can be drawn using the elements i) - iii) mentioned above) to render the model finite. 

Moreover, the counterterm corresponding to this block is sufficient to regularize all higher 
moments of fields u and m', since any diagram of perturbation series is expressed as a convolution 
of equivalent blocks. Diagrams contain no additional superficial divergences [TB]. We recall 
that all loops which could arise in the diagram expansions are created by the single local vortex 
with any number of m'— legs incident at it. The counterterm corresponding to the elementary 
divergent block is constant and local in configuration space, i.e. oc bit — to)'^(x — xq). In 
the action functional ([H]), the same local term appears, so that the model is renormalized 
multiplicatively, and only the renormalization constant Z is required. To keep the interaction 
coupling constant dimensionless, one must introduce a mass parameter [i. 

The bare parameters are related to the renormalized parameters by 



The auxiliary fields and the Green function are related to their renormalized analogs by: 
The only renormalization constant required in the model is 



1=1 

where the amplitudes ci are defined to be precisely those needed to subtract the poles in the 
corresponding Feynman integrals. The derivation of the RG equation is carried through in 



r]o = w Z' 





(12) 
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a standard way since the curvature coupling parameter ^ is not renormalized. We apply the 
standard RG differential operator nd^ (for fixed uq, rjo, ^o) to Gr to obtain 



d di] d 



^ 1 7 



(13) 



It follows from the relation (ITT]) that for the RG functions (the coefficients in the RG equation) 
one has 

7(r/) =/i^lnZ, /5(ry)=r/(-2£-(a-l)7(77)). (14) 



In the forthcoming section, we discuss the existence of a fixed point of RG transformation stable 
with respect to the long time large scale asymptotic behavior, that is a solution rj^ of equation 
P{ri) = such that I3'{rj^) > in the physical range of parameter rj > 

In the equation of critical scaling for the curved space, one has to take into account the 
new mass parameter w in addition to ?7, ^, and /i. It had been introduced in [38] to quantify 
the dilatation w'^gij of the background metric Qij. The derivation of the scaling equation is 
standard: the canonical scale invariance of the renormalized Green function Gr{C,, /x, z/, rj) with 
respect to dilatations of all variables is expressed by the equations 



d w 



-2z/- X- d 

aw Of ox 



(15) 



dt dv 



In the IR fixed point of RG-transformation, one excludes the differential operations fidfj, and 
udi, using the RG equation (IT^ to obtain the equation of critical scaling. 



-zxj- h rf - 7 



GR{^,fi,i^,v) = 0, 



(16) 



where vu is now interpreted as the metric scaling variable. If a stable point of RG transforma- 
tion relevant to the long range asymptotical behavior 77^, exists in the model, the value of the 
anomalous dimension 7 at the fixed point is found exactly, owing to the relation between f3 and 



7, 



7* 



2e 



1 - 



5x 



6^ — d. 



" = "lo 



a — 1 a — 1 \ d 

In agreement to dimensional considerations, the renormalized Green function has the form 

.2 1 



(17) 



GR{t,r) = {ut)-''/\ 



— .wvt 



x — x 



where x is a scaling function of the dimensionless variables. The dependence on r] is not 
displayed explicitly, because the derivatives with respect to this parameter do not enter into 
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the scahng equation. It follows from (ITBI) that at the fixed point x satisfies the equation 



' _^ _ 7* 
^dz 2 



x{s,y,z) = 0, 



its general solution is xi^iVi ^) = s'^*^'^ijj{s,y) where tp is an arbitrary function of the first and 
second arguments. For the Green function ffTOl) we then obtain 

G{t, r) ~ t-'^/2-7,/2^ 1^^^ ^^^^ _ t-^^/^^jj 1^^, caz/t^ , 5, > d, (18) 

where the form of the scaling function ip is not determined by the equation flTBl) . Although the 
value of 7* in ( flTl) and the solution ( fTSl) have been obtained without calculation of the constant 
Z, such a calculation is necessary to check the existence and stability of the fixed point. Within 
the e— expansion these facts can be verified already in the one-loop order. We perform it in the 
forthcoming section. 

We conclude the section with a remark on the arguments of scaling function ip for the planar 
3-regular graphs of order 2N with the standard orientation (i.e. the 2A^-honeycomb). In the 
Sec. 2, we have mentioned that they are equivalent to the sphere of radius p. The surface area 
of the sphere equals to 2'kN = 47rp^, and therefore p = a/ N/2. The relevant Gilkey coefficients 
and the asymptotic behavior of the Feynman propagator as t — > are given at the end of 
Appendix A. One can see that the corrections to the standard diffusion kernel risen by the 
space curvature can be naturally interpreted as the corrections due to the finite size of the 
3-regular subgraph. Then, the scaling function ip is 

and it can be calculated by means of diagram expansion. In the thermodynamic limit, N ^ oo, 
(when the the graph F is large and regular) these corrections vanish and the scaling function 
depends only on one argument r^/ ut relevant to the fiat space. The result on the critical scaling 
( IT8|) is still valid in the thermodynamic limit. 



ip = ip 



7 A fixed point of RG transformation and heat-kernel expansion 

The main technical difficulty concerning diagram computations in curved spacetime is the 
lack of a global momentum-frequency representation. The momentum space is associated to 
each point x' by the Fourier transformation 



f{x,x') = 
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where ky = g^^kiUj which does not help much in calculations. In the flat metric, the calculation 
of the residue ci in (fT2l) . in the one- loop order, can be performed in the standard way using the 
momentum-frequency representation 



c,U, = ^a-V("-), Z = l + |. (19) 
It corresponds to the fixed point 

which is positive and stable for small e > and a > 1. In the curved space the new singularities 
would arise at small scales (as x x') even for the metrics of constant curvature. They do not 
affect the renormalization group procedure (at least at the one-loop order) at large scales, but 
breaks down the perturbation theory for the even space dimensions. 

The necessary regularization of effective action can be performed by using the heat-kernel 
expansion which does not depend on the spacetime topology (see, for instance, [3S|). We 
consider the operator Q2 with the kernel 

where Sn is the renormalized action functional of the model we consider. The one-loop amend- 
ment to the renormalized action is given by 

Li = -^Trln(g2A-i) (21) 



where the propagator matrix is given by 

A = 



A 

A,,/,, 



UW 



The counterterm corresponding to the only divergent diagram block is linear in the auxiliary 
field u'. Variation of Li with respect to the auxiliary field 

6^>L^ = ^Tr {AQ,'6^iQ2A-')) , S^iQ^A-') = -a{a - l)7]uu'''' A-\ 



can be written as 



5u'Li = --5u' |Tr J — exp(-r ■ Q2A"^)| 



Let us suppose that A(t, x, y,(52A ^) is the kernel associated with the operator exp(— r 
Q2A~^). Then, the one-loop contribution fl2T|) reads as 



2 r roo 

^1 = -- J dv^ J — A(r, x, X, QsA-^) 



(22) 
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Since the divergences are purely logarithmic, one can neglect the inhomogeneity of the scalar 
field M in X while calculating the relevant counterterm. Being treated as a constant, the scalar 
u renormalizes the scalar curvature, 

so that one can define the modified Feynman propagator, Auu' ^'uu'- 

Then, it appears that the kernel A(r, x, x, (A')~^) meets the "heat equation" 

[9. + (A')-^]A(r,x,x,(A')-^) = 

and therefore allows for the asymptotic expansion as r 0+, in the limit x x', 

oo 

A(r, X, X, ~ (47rr)-'^/2 ^ r^^E^^, A'), (23) 

m=0 

in which A') are the relevant Gilkey coefficients |10]-[l3]. In the case of a constant 

curvature metric, the first coefficients can be found in the Appendix A fl27j) . Substituting the 
expansion (1231) into the integral representation (1221) . one can perform the integrations with 
respect to the parameter r with a cut off tq. As a result, one arrives at the power series, 

00 m-d/2 

-(4?)^£;?ri^75^.»(x.A'). (24) 

being singular for integer m such that d = 2m. It is obvious that the breakdowns of perturbation 
theory at small scales would happen for the even space dimensions. Extracting the pole part 
of Li one finds that 

where -Erf/2 are the relevant Gilkey coefficient provided that d is even. It is important to mention 
that in (1251) the singularity coincides with that of (ITQll . The residue in ( l25ll has the "good" 
signature and does not break the stability of the fixed point of RG transformation. 



P.P.[Li] = -jj^ I dv,Ea/2{x) = -jj^ I dv^Ea/2{x) (25) 



8 Discussion and Conclusion 

In the present paper, we have studied the transport through the large complex networks 
containing regular subgraphs. We have considered such networks as generalized trees in which 
two types of nodes are allowed: simple nodes and supernodes. We have supposed that the 
supernodes are either the subgraphs with many polygons or A;-regular subgraphs. In particular. 
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we have discussed the case of 3-regular subgraphs which can be treated as the complete Riemann 
curved surfaces characterized by finite areas and can be compactified. 

The diffusion process taking place on such a complex network is considered as a generalized 
Brownian motion with arbitrary boundary conditions. Its random dynamics is then described 
by a functional integral. We have studied the long time large scale asymptotic behavior of 
the Green function for the nonlinear diffusion equation defined on the complex network sup- 
plied with an integrable initial condition. The transport trough the complex network is of a 
strongly nonlinear nature being affected by both the varying effective space dimension between 
supernodes and the space curvature within a supernode. 

In the regular graphs endowed with the standard orientation of edges (the cyclic ordering 
of edges to be the same for all nodes), the effect of curvature within the supernodes reveals 
itself by the finite size corrections to the scaling behavior. In the space of even dimensions , 
the curvature results in the additional singularities of Green function at small scales and calls 
for the special regularization. We have to stress that in contrast to the case of gravitational 
field, we are not restricted by the equivalence principle while considering models of complex 
networks. The main technical difficulty of computations in the curved space metrics is the lack 
of the global frequency- momentum representation. We have used the heat-kernel expansion 
which does not depend on the space time topology to regularize the effective action at small 
scales in the one-loop order. 

Our approach to the complex networks containing large well-structured regular subgraphs 
could be used in studies of spreading of viruses [52] or social epidemics [53], traffic properties 
and the modelling of various ecological webs. 
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Appendix A. Feynman propagator in curved space 

We present the explicit solution of the linearized diffusion problem, in the d-dimensional 
curved space metric, for G(x, x', t) satisfying limx'-*o G'(x, x', t) = G(x, t) and limt_+o G'(x, x', t) = 
5(x — x'). The details can be found in [IB] for field theories and in [U] (see also references 
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therein) for classical models. 
The general solution is 

A„„,(x,x',t) = g(t) "^P^^;f^ xp (-%^) AV2(x,x')fi(x,x',t) (26) 

where 6{t) is the Heaviside function, cr(x, x') equals to half the square of the geodesic distance 
between x and x', and A(x, x') is the van Vleck determinant, 

det(pa(x^ 

which reduces to unity in flat space, R is the scalar curvature. The function f2(x, x',it) allows 
for the following series expansion in the limit x' x: 

oo 

lim ^](x, x', t) = 1 + V t^Ei{yi) 

1=1 

valid in the limit t — > where -E'z(x) are known in the literature as Gilkey coefficients |40]-|43j. 
For the diffusion equation, in the flat metric, the only coefficient which contributes is Eq and 
we recover the well known standard diffusion kernel. 

The planar 3-regular graph of order 2N with the standard orientation (the cyclic ordering 
of edges is taken the same for each node), a 2A^-honeycomb, corresponds to a sphere of radius 
p = \JN/2. The Ricci scalar curvature is then 

and the Gaussian curvature equals to k = 2/N. The Gilkey coefficients reduces to 

2 

Eo = l, E, = —, 

^^-isiv^' ^'-3iEW^- ^^^^ 

Then, in the limit x' ^ x and t 0, the Feynman propagator exhibits the following dependence 
on the size of 3-regular subgraph: 

^^P / 1 f2uot\ 1 /2i/otV 4 (2vQt^^ 
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Figure 1: In the model, we consider the complex networks as the generalized trees with two types of nodes: 
the simple nodes (probably of low connectivity) and the supernodes which are either the subgraphs with many 
triangles or the fc-regular subgraphs. 
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Figure 2; Three diagrams of the Feynman diagram expansion for the Green function of the nonlinear diffusion 
equation. 



23 



